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Abstract. 

We reevaluate a quark model prediction using the new QCD evolution function calcu- 
lated to the 3 loop order and conclude that this model compares favorably with the new 
experimental results. 
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From the semi-leptonic decays of baryons the matrix element of the flavour changing 
parts of the axial currens are fairly well known |2[ . The first EMC measurements of the 
spin dependent proton structure function g^{x,Q'^) when integrated, gave information of 
the matrix element of the flavour singlet part of the axial current between proton states. 
The first experimental results Q 

rp(Q2) = f dxgP{x, Q2) = 0.126 ± 0.012(stat) ± 0.026(syst), (1) 



implied that the matrix element of the flavour singlet part of the axial current between 
proton states was compatible with zero within the error bars. Despite the large experimen- 
tal errors, some models of the baryon structure would give results outside three standard 
deviations of the measured flavor singlet matrix element value. Since this was an unex- 
pected result, it generated a lot of discussion, see a short review in Ref. and forced 
us to reconsider our present model understanding of the baryonic structure. However, the 
Skyrmion model of Ref. [|| and the chiral bag model Q was immediately recognized to give 
results that was within one standard deviation of the measured integral over g^{x,Q'^). 

In the last two years there has been a lot of new important experimental results com- 
ing from deep inelastic polarized lepton scattering both on polarized proton, deuteron and 
helium targets |^, |l^ . It is of interest to see how well these early model calculations 
compare with the new experiments when we also take into account recent progress in the 
calculation of QCD corrections. The spin structure function g\{x,Q'^) has now been mea- 
sured at quite small values of x, and a smooth extrapolation to x=0 Q gives a new value for 
the integrated proton spin structure function Tp{Q'^) at = 10 GeV'^ /c^ Q. 

rp(Q2) = 0.136 ±0.011 ±0.011 (3) 

This value is larger but compatible with the earlier result, eq.(|l]). To compare this new 
value with model prediction, requires a good knowledge of QCD radiative corrections and 
in particular the SU(3) flavor singlet QCD evolution function, which, it was argued, may 
strongly affect the quark model predictions |l^]. Very recently this evolution function has 
been evaluated at the 3-loop level [^] , and confirmed by Larin [0] , who gives references to 
earlier works and who also refers to the order QCD radiative corrections calculations. 

In the following we shall argue as if nature is SU(3) flavour symmetric, not because we 
believe that to be an extremely good approximation, but because frequently this assumption 
is used in extracting matrix elements of currents with a definite flavour from the experimen- 
tal semi-leptonic decay data |[l|, 0. Two main points are responsible for the present model 
understanding of the baryon structure Tp{Q'^). The first is the role of the axial anomaly 
which we shall discuss in the chiral SU(3) limit. The second point concerns the baryon wave 
function itself and what we can learn from the other baryon spin observables. 

To discuss the first point we present our model picture of the baryon structure. We 
shall assume that an effective Lorentz-scalar force confines the quarks to a finite region in 

^See a recent discussion on this point by Close and Roberts pill 

^ The SM CoUaboration have also evalated a "world average"0| at = lOGeV^ of 

Tp(lQGeV'^) = 0.142 ± 0.008 ± 0.011) (2) 
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space (which we shall take to be a spherical cavity in the rest system of the baryon) where 
they interact through the exchange of gluons. As the valence quarks are excluded from the 
QCD vacuum outside the cavity, they will polarize the region just outside the cavity in qq 
pairs in such a manner that the axial current carried by the valence quarks extends into the 
QCD vacuum surrounding the cavity by qq pairs carrying the quantum numbers = 0~ , 
i.e., the pairs are the effective Goldstone bosons of the QCD vacuum. For chiral symmetry 
to be a good symmetry, it is necessary (but not sufficent) that the effective Goldstone fields 
in the QCD vacuum couple to the axial currents of the quarks in the cavity such that the 
axial current is continuous in all space. 

In the chiral SU(3) limit the valence quarks (u,d,s) in the cavity are massless. The 
pseudoscalar Goldstone bosons surrounding the quark confinement region, i.e., the pion, 
the rjs and the kaon, are also massless. (This naturally imply that the qq content of the 
proton is mainly a large distance phenomenon.) On the other hand, the flavour singlet 
pseudoscalar, the rj', is massiv in the chiral limit due to the U(l) anomaly. It is this 
consequence of the anomaly that we shall explicitly use to illustrate how a polarised ss 
cloud could exist around the "bare" nucleon (the nucleon's quark core), even if the "bare" 
nucleon contains no strange quarks. To illustrate this we use the following simplistic picture: 
Suppose for the sake of argument that we have a single confined u quark. Since we are in 
the chiral limit the axial current, which is carried by the u-quark inside the confinement 
cavity, has to be continued into the vacuum outside the cavity by massless pseudoscalar 
mesons, which would carry the flavour content of a uu state. This uu state can be written 
in terms of states that have well defined transformation properties under SU(3) flavour at 
the confinement "surface" where it changes into pseudoscalar mesons: 



If all pseudoscalars were massless, the ss content of rjs and r]' would cancel trivially, as is 
implied by this equality. However, the three pseudoscalar mesons have different masses, 
and therefore their Yukawa ranges beyond the cavity surface will be different. For example, 
even in the chiral limit the heavy rj' (heavy due to the anomaly) and the massless T/g 
would generate a strangeness content in the proton since the strength of the heavy rj' field 
would decay exponentially outside the confinement "surface" . Therefore, the "bare" nucleon 
containing only nonstrange quarks has a cloud of ss like states even in the chiral SU(3) limit 
as a consequence of the U(l) anomaly. 

The second important point in these considerations is as follows: If one wants to 
compute matrix elements of the flavour singlet current, there is no reason to have confidence 
in the results if the model is not giving correct results for the matrix elements of the 
flavour octet currents which are responsible for weak semi-leptonic baryon decays. From the 
experience with the MIT bag model |]T5| we know that relativistic quarks reduce the neutron 
axial charge gA from 5/3 towards a more realistic value. In addition, it is necessary to go 
beyond the usual zero order SU(6)-like baryon wave functions of an additive quark model. 
An additive quark model gives (in the chiral limit) for hyperon P- decays an amplitude 
ratio, F/D = 2/3 like the non-relativistic quark model. The fallacy of the additive quark 
models becomes apparent for the baryon magnetic moments where these quark models give 
the wrong inequality, i.e., for the measured magnetic moments |//(H~)| > |^(A)| |^]. These 
"spin" problems of F and D and of the baryon magnetic moments are easily explained by 
invoking the same effective color magnetic spin-spin interaction that give the " bare" nucleon 




(4) 
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- "bare" A mass splitting 17, T^. This effective interaction introduces dynamic spin-spin 
correlations in the baryon wave functions that not only resolves the F/D problem, F/D ~ 
0.57, i.e. < 2/3, but also easily explains why > |^(A)|. This spin-spin correlation 

contributes to the neutron's gA and is important for the axial flavor-singlet charge 



As should be evident from the above discussion, to make a sound model prediction 
for Tp{Q'^) and r„(Q^), it is necessary that the quark model reproduces the value for 
gA{= 1.2573 lb 0.0028) of the Bjorken sum rule as well as other baryon spin observables. 
Then one can discuss the model prediction for the flavor singlet part of the integrated proton 
spin structure function and the so-called "spin content of the proton". To be spesific we 



first define the relevant flavor SU(3) amplitudes P, 2C]. For Q — > oo 



rp(Q2 ^ oo) =< p T I E \^i^Hi^'^i\p T> = \[h + h + h] (5) 

i 

where qi is the charge of the i'th quark. The matrix element of the axial current operator 
between proton states on the r.h.s. of eq.(|5|) is evaluated at = 0. To make contact with 
the experimental result of eq.(P) the axial current operator on the rigth-hand-side of eq.(^) 
should be renormalized at the point = 10 GeV"^ In eq.(^) the SU(3) flavor singlet 
amplitude is 

= '\/|<pT |%573Ao^|pT> (6) 



and the two flavor octet amplitudes are 

1 



° 2^3 

and 



<pT |^7573A8'I'|pT> (7) 



^3 = ^ <pT I%573A3^'|pT> (8) 

where A are the usual Gell-Mann matrices. In terms of the standard SU(3) amplitudes F 
and D for the baryon semi-leptonic decays, the two flavor octet amplitudes are 

/3 = \{F + D) and h = \{F - D/3). 

The correspondence between the /j's and the commonly used Ag's are: 

Au = i/o + Is + h, 
Ad= I/0 + /8-/3, 
As = i/o - 2/8 , 

and the "proton spin content" is: 

S = Au + Ad + As = |/o 

Since the two amplitudes I3 and /g are matrix elements of conserved currents in the 
chiral SU(3) limit, their values are independent of the renormalization point. Therefore, 
models used to calculate F and D {I3 and /g) should have conserved axial octet currents 
in the chiral limit in order to be reliably used at the scale of the SMC experiment. This 
is a key requirement in any model prediction of the " proton spin content" , a requirement 
which is not satisfied by the non-relativistic quark model, for example. A concrete example 
of a model satisfying the above condition is the chiral bag model, an extention of the 
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MIT model. Calculations within this model are quite straigtforward, and give for many 
observables excellent insight and estimates of the relevant physics involved. This was (and 
is) for us a valid reason for our use of this model. One of the serious drawback of this 
and similar models is the sharp confinement "wall" (like a square well potential) and the 
associated translational invariance problem. 

Before presenting our quark model results it is necessary to outline the radiative QCD 
corrections to this sum rule and the possible renormalization of the flavor-singlet axial 
current. The proton sum rule including the QCD corrections at the scale have recently 
been calculated at the 3-loop level of perturbative QCD (for references see and the 
later next order guestimates of Kataev pi]]) and reads, 

r,(Q^) = 1[1 - ^ - 3,5833 f 2(«!lV - 20.2153 f " 1™ 

6 TT \vr/ \vr/ 




+ - - 1.0959 f 2<5!)y _ 3.7 f 515!))'] (9) 

As discussed, this expression assumes that the two flavour octet amplitudes, and Ig, are 
independent of the scale. On the other hand the scale dependence of the flavour singlet 
amplitude, /o(/U^), is made explicit in this equation. The effect of the anomalous dimen- 
sions of the flavor-singlet matrix elements are taken care of by the renormalization group 
(exponent) evolution function y(a(/i^), a((5^)), which very recently has been evaluated at 
the 3-loop level (in the lowest twist approximation), e.g. [|l4|, and reads 

lnY{a{,ix\a{Q')) = - \ da'- / - (10) 



where 7^^) = (1 - n) + 11 /?o = H - 2 n//3. Pi = 102 - 38 n//3, 

and P2 = - ^xl^'^/ + ^™/- The integrand is an increasing function of a' in the region 

where the formula is valid. 

To define the renormalization group-invariant (i.e., convention independent) nucleon 
matrix element of the singlet axial current |jl^, ^ , we will write the evolution function as 
a product Y{a{fi'^),a{Q^)) = Y{0,a{Q'^))Y{a{n^),0). Then 

^inv = exp ( -— "^"'^ J ^(-"'^ = y(a(A.2),0)S(/.2), (11) 

and since a(Q^) for = 10 GeV'^ jc^ is small, we can use eq.(p^ and write the expansion 
for nj = 3 

y(0,a(Q2))~ 1 + 0.66667^^ + 1.213 1^^1 +••• (12) 

TT \ vr 
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It then follows that the QCD corrections for the integrated spin-structure function are 
contained in QCD coefficients that are convention independent |14, 21 1: 




r,«3=, = i|l - - 3.5833 (^) ' - 20.2153 (^) ' - 130 ^ ' 



In the following we shall discuss the influence of the flavor singlet evolution function 
on the sum rule in eq.(|^) and how the flavor singlet amplitude lo^modeh eq.(|6|) (or Smodei) 
calculated in the usual spectroscopic quark models which include quark conflnement, should 
be included in the two expressions, eqs.(|9|) and (|l3|). Two extreme viewpoints regarding 



the influence of the flavor-singlet evolution function Y {a{^'^) , a{Q^)) |T^, ^] and |23] were 
discussed, and still is a matter of debate. One viewpoint argues that it is not inconceivable 
that non-perturbative effects could cause an evolution down to the very low spectroscopic 
scale //^ such that Y{a{ii^),a{Q'^)) is a very small number, a number which in itself could 
explain the EMC result ]T2| , |22[ . This scenario implies that the conflnement quark model 
calculation predicts the matrix element of the flavour singlet axial current /o(/^^) in eq.(|9|). 
Furthermore, such an argument would lead to practically no constraint on a quark model 
calculation of /q. In one of our previous works we adopted this viewpoint, see e.g. 
Refs.[12, 24|, and discussed the effects the known behaviour of the renormalization group 



exponent could have on the "proton spin content". This contrasts with our original atti- 
tude : using models with conflned quarks to compute the matrix elements appropriate for 
the Bjorken limit, we compute physical observable quantities, quantities which should be 
independent of the renormalization point and therefore should be identifled with the renor- 



malization group invariant quantities |23] or Yjijiy of eq. (|ll| ) . This means only perturbative 
QCD corrections outlined in eq.(|l3|) should be included in the sum rule. The reason being 
that any realistic quark confinement model (including chiral bag models) has built-in long 
distance effects, meaning soft gluonic corrections are to some extent included in the quark 
wave functions of these models. In addition, the chiral bag model also includes the "very 
soft" effective mesonic mediators of the longest range nuclear forces. These effective qq 0~ 
degrees of freedom are important for understanding properties of baryons. This is of course 
very different from the parton model approach where in the first appoximation the valence 
quarks are free non- interacting quarks with no gluonic content. 

After these discussions we can present some numerical estimates of what is expected for 
the term y(a(/x^), a((5^)) Io{^J^) of eq.(P) and compare with the ^modei value calculated in 
the chiral bag model where the pseudoscalar mesons are excluded from the bag. We use the 



values for F and D of the chiral bag model calculated in the chiral limit 1 2C ] including the 
two-body magnetic gluon exchange spin-spin correction: F = 0.455 and D = 0.795 {F/D 
= 0.57). This means that in eq.(P) Is + Ig = F + L>/3 = 0.72. As values for a{Q'^) we use 



the ones corresponding to A=240 MeV [P5| that runs within one standard deviation of data 
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from the Zg-mass to the r-mass. Then we find using the experimental values of eq.(^): ^ 
y(a(^^), a(Q2))/o(/i^) = 0.194 ± 0.073 ± 0.073. (14) 

This result should be compared with the direct calculation of the chiral bag, which gave 
Io,modei = 0-22 corresponding to a "spin content" of SmodeZ = 0.33 when the pseudoscalar 
mesons were excluded from the interior of the bag. ^ Early in the spin-crises period |^] , this 
calculated value, T^modei = 0.33, was considered to be too large, which certainly, with the new 
experimental results 0] and the refined calculations of the radiative corrections, no longer 
is the case. The latest Elhs and Karhner [||] analysis gives S(Q2 = IQGeV"^) = 0.31 ±0.07 
(when higher twist effects are ignored). Since lo^modei is very close the the present value of 
eq. (|l|), there is in this model little room for the scenario with Y being very small. However, 
as discussed above, since the quark wave function of the chiral bag model contains soft-gluon 
effects, we should and shall identify Smode« = 0.33 with of eq.(|l3|) when we calculate 
Tp{Q'^). The corresponding result for r„((5^) is obtained from the above by substituting 
Is ~ I3 = -0.53 for /g + /s = 0.72 in the same equation (|l^). In figures 1 and 2 we show 
both the model results using eq.(p^ and the experimental data as a function of a((5^). We 
have not included any corrections from higher twists that would influence the results at low 
Q^, the reason being that there is doubt both about the size as well as the sign of these 
corrections [^]. As we can see, the model is giving results that are very satisfactory when 
compared with the experiments on both the proton and the neutron targets. 

In this model the so-called proton spin content has many components. One third comes 
from the spin of the three valence quarks, the remainder originates from the orbital angular 
momentum (of the relativistic quarks) and from the gluonic exchanges between the valence 



quarks which can be regarded as confined qq pairs |15]. Note that the U(l) anomaly is 
essential in reconciling the correct values of F and D with a value of AS < 0.57. In terms 
of the Ag's the model gives An = 0.83, Ad = -0.42 and As = -0.08 as 00. 

We have also noted that a recent calculation by Narison, Shore and Veneziano based 
on QCD spectral sum rules give results ||2^ that are very similar to the ones we have 
presented. Our results were calculated in the chiral limit for the three quarks where SU{3) / 
by definition is a good symmetry group. When quarks and pseudoscalar mesons are given 
masses (m„ = rrid ~ 10 MeV and mg « 200 MeV), SU{3)f is broken, then what are the 
changes? The neutral pion is fairly light but both rj and rj' become heavy. Due to their 
heavy masses, the influence of both of the two pseudoscalar ss carrying mesons are then 
suppressed. This leaves almost no hidden strangeness around the nucleons, as shown in an 
earlier calculation for nucleons Q. However, the results for rp^„,(oo) hardly change at all!, 
meaning the chiral bag model with broken SU{3)f symmetry and almost no ss content can 
describe data as well as the SU (3) / symmetric chiral bag model above. From these results 
it follows that the experimental data of today do not necessarily imply that As 7^ but 
but could just as well mean that SU{3) / symmetry is broken. This specific example should 
underline the fact that the extraction of As from the data on the basis of SU{3) / symmetry 
should be taken with some caution. 



^ If we use the "word average" value 

Y{a{fi^),a{Q^))Io{i/) = 0.213 ± 0.053 ± 0.073. 

* In a version of the model where the pseudoscalar field is continuous through the bag surface (cloudy 
bag), lo^cioudy = 0.38 corresponding to T.cioudy = 0.57 Q, the same value as in the Ellis and Jaffe sum rule. 
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When a model is able to reproduce pre-EMC data on static properties of the baryons 
such as the axial charge, coupling strength of pions to the baryons, magnetic moments etc., 
and also give reasonable fits to the SMC data in a new domain when the effect of the U(l) 
anomaly on the rj' mass is taken into account, there is a coherence that give us confidence 
that we have some understanding of the rudimentary aspects of the nucleon structure. 

Correspondence with E. Reya are gratefully acknowledged. This work is supported in 
part by NSF grant no. PHYS-9310124, a NATO travel grant and a grant from Norges 
Forskningsrad. 
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1 Figure Caption 

Figure 1: Comparison of the chiral bag model calculation, eq.(^), with experimental 
data of Tp{Q'^). Systematic errors are heavy vertical lines, statistical errors are 
thin vertical lines. Data are plotted for reasonable values of a{Q'^). 

Figure 2: Same as Figure 1 but for r„(Q^). 
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